Monolayer honeycomb structures of group IV elements and III-V binary compounds 
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Using first-principles plane wave calculations, we investigate two dimensional honeycomb structure 
of Group IV elements and their binary compounds, as well as the compounds of Group III-V 
elements. Based on structure optimization and phonon mode calculations, we determine that 22 
different honeycomb materials are stable and correspond to local minima on the Born-Oppenheimer 
surface. We also find that all the binary compounds containing one of the first row elements, B, 
C or N have planar stable structures. On the other hand, in the honeycomb structures of Si, Ge 
and other binary compounds the alternating atoms of hexagons are buckled, since the stability is 
maintained by puckering. For those honeycomb materials which were found stable, we calculated 
optimized structures, cohesive energies, phonon modes, electronic band structures, effective cation 
and anion charges, and some elastic constants. The band gaps calculated within Density Functional 
Theory using Local Density Approximation are corrected by GWo method. Si and Ge in honeycomb 
structure are semimetal and have linear band crossing at the Fermi level which attributes massless 
Fermion character to charge carriers as in graphene. However, all binary compounds are found to be 
semiconductor with band gaps depending on the constituent atoms. We present a method to reveal 
elastic constants of 2D honeycomb structures from the strain energy and calculate the Poisson's 
ratio as well as in-plane stiffness values. Preliminary results show that the nearly lattice matched 
heterostructures of these compounds can offer new alternatives for nanoscale electronic devices. 
Similar to those of the three-dimensional Group IV and Group III-V compound semiconductors, one 
deduces interesting correlations among the calculated properties of present honeycomb structures. 

PACS numbers: 73.22.-f, 61.48.De, 63.22.-m, 62.23.Kn 



I. INTRODUCTION 

Last two decades, nanoscience and emerging nan- 
otechnologies have been dominated by honeycomb struc- 
tured carbon based materials in different dimension- 
ality, such as fullerenes, single and multi walled car- 
bon nanotubes, graphene and its ribbons. In particu- 
lar, graphene, a two dimensional (2D) honeycomb struc- 
ture of carbon, has been an active field of research. 1 
Because of unique symmetry, electron and hole bands 
of graphene show linear band crossing at the Fermi 
level 2 - resulting in a massless Dirac Fermion like be- 
havior of charge carriers. As a result, Klein para- 
dox, an interesting result of quantum electrodynamics 
was expected to be observed in graphene More- 
over, it was shown that half-integer quantization of Hall 
conductance^^ can be observed in graphene. Unusual 
electronic and magnetic properties of graphene, such as 
high carrier mobility, ambipolar effect, have promised va- 
riety of applications. In addition to some early works 
on crystalline order in planar structures'* 1 ^ possibil- 
ity of very large one-atom-thick two dimensional (2D) 
crystals with intrinsic ripples is reported theoretically-ii 
and experimentally^ 2 - Not only extended 2D graphene 
sheets but also quasi-one-dimensional graphene ribbons 
with armchair or zigzag edges have shown unusual 
electroniC ) 13 i 14 ' 15 i 16 i 17 i 18 magneti o 19 i 20 ' 21 i 22 and quantum 
transport properties. 23 ' 24 ^^ 

All these experimental and theoretical studies on 
graphene created significant interest in one-atom-thick 



honeycomb lattices of other Group IV elements and com- 
pounds of III-V and II- VI Group elements. Recently, the 
boron-nitride (BN) honeycomb sheet was reported as a 
stable ionic monolayer.— BN has the same planar struc- 
ture as graphene with a nearest neighbor distance of 1.45 
A. However, its ionic character causes a gap opening at 
the K-point. Thus, instead of being a semimetal, BN 
honeycomb structure is a wide band gap insulator with 
an energy gap of 4.64 eV. Soon after its synthesis, several 
studies on nanosheets2£ and nanoribbon a 29 ' 30 i 31 i 32 of BN 
have been reported. 

Hexagonal monolayer of Zinc Oxide (ZnO) is a II- 
VI metal-oxide analogue of graphene and BN. Pre- 
viously, works on nanostructures of ZnO such as 
nanosheets' 3 - nanobelts^ nanotubes^ nanowires 3 ^ and 
nanoribbons^ were reported and recently the synthesis 
of ZnO bilayer honeycomb structure was also achieved. 38 
In contrast to graphene, ZnO nanoribbons have fer- 
romagnetic order in their ground state due to elec- 
tronic states at the zigzag edges dominated by oxygen 
atomsi 39 ' 40 

Two dimensional SiC honeycomb sheet is another 
Group IV binary compound displaying interesting prop- 
erties. While the infinite periodic 2D form of SiC is a 
semiconductor with 2.55 eV band gap, and its zigzag 
nanoribbons are magnetic metals, the armchair ribbons 
are nonmagnetic semiconductors. 41 Half-metallicity is 
also predicted for narrow SiC zigzag nanoribbons without 
any chemical decoration or applied external field. Fur- 
thermore, functionalization of SiC single sheets upon for- 
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FIG. 1: Top and side views for two dimensional (a) pla- 
nar (PL) and (b) low-buckled (LB) (or puckered) honeycomb 
structure. In the PL structure atoms are located on the same 
plane. In the LB structure the alternating atoms are located 
in two different parallel planes. The buckling A is the dis- 
tances between these two planes. Bravais lattice vectors for 
both structure are given with |oi| = |ct2 1 = a. The unitcell is 
delineated and shaded. 

mation of various types of vacancies and adatom decora- 
tion was also predicted.— 

Very recently, we have reported that among Group IV 
elements, not only C but also Si and Ge can form stable 
honeycomb structures.— It is found that for Si and Ge 
planar (PL) geometry is not the lowest energy configura- 
tion and it is not stable. Alternatively, it was shown that 
a low-buckled (LB) (or puckered) geometry corresponds 
to a stable local minimum on the Born-Oppcnhcimer sur- 
face. Buckled honeycomb structure of Si was pointed out 
even in some earlier studies ^ 44 ' 45 Surprisingly, in spite 
of their puckered geometrical structure, Si and Ge mono- 
layers have electronic band structures which are similar 
to graphene. As a result, linear crossing of tt- and ir*- 
bands at K- and K'-points of the hexagonal BZ attributes 
a massless Dirac fermion character to the charge carri- 
ers. Quasi ID honeycomb structures, namely nanorib- 
bons of Si and Ge, also show interesting electronic and 
magnetic properties depending on their width and ori- 
entation. Successful realization of single crystal silicon 
monolayer structure o 46 ' 47 through chemical exfoliation 
shows that 2D silicon monolayers with their low resistiv- 
ity and extremely thin structures can be quite promising 
for nanoelectronics. 

Motivated by the recent experimental developments 
and theoretical investigations on 2D monolayer honey- 
comb structures, in this paper we carried out a system- 
atic study of similar structures of Group IV elements 
and III-V binary compounds based on first principles 
calculations within Density Functional Theory (DFT). 
Our objective is to reveal whether monolayer honeycomb 
structures can be found as a local minimum on the Born- 
Oppcnheimer surface. We hope that the predictions of 
our work will guide further experimental studies towards 
the synthesis of new materials having honeycomb struc- 



ture. The present work, which considers a total of 26 
elemental and binary compounds in 2D honeycomb struc- 
ture and reveals whether they are stable, is an extension 
to our preliminary work on Si and Ge puckered honey- 
comb structures^ Based on extensive analysis of stabil- 
ity, 22 different materials out of 26 are found to be stable 
in a local minimum on the Born-Oppenheimcr surface ei- 
ther in finite size or in infinite periodic form. We hope 
that interesting properties predicted by this study will 
promote efforts towards synthesizing new materials and 
heterostructures, which will constitute a one-dimensional 
analogue of 3D family of tetrahedrally coordinated semi- 
conductors. 

The organization of this paper is as follows; In Sec. II 
the methods together with parameters used in our calcu- 
lations are outlined. In Sec. Ill, we determine the atomic 
structure and related lattice constants of the honeycomb 
structures via total energy minimization. We also discuss 
how the stability of the structure is maintained through 
puckering. In the same section, we present our results 
regarding the calculation of phonon modes and our anal- 
ysis of stability based on these results. The mechani- 
cal properties of these structures are investigated in Sec. 
IV. We discuss the electronic band structure of various 
stable materials calculated within DFT in Sec. V. The 
underestimated band gaps are corrected by using GWo 
calculations. As a proof of concept for a possible fu- 
ture application of these materials we consider semicon- 
ductor superlattices formed from the periodically repeat- 
ing pseudomorphic heterostructure in Sec VI. We showed 
that the superlattices have an electronic structure differ- 
ent from those of constituent materials and behave as 
multiple quantum well structures with confined states. 
In Sec. VII, our conclusions are presented. 



II. METHODS 

We have performed self-consistent field, first-principles 
plane-wave calculation s 48 ' 49 within DFT for total en- 
ergy and electronic structure calculations. Projec- 
tor augmented wave (PAW) potentials^ and exchange- 
correlation potentials approximated by local density 
approximation 51 (LDA) are used in our calculations. In 
the self-consistent field potential and total energy calcu- 
lations a set of (25x25x1) k-point sampling is used for 
Brillouin Zone (BZ) integration in k-space. Kinetic en- 
ergy cutoff ?L 2 |k + G| 2 /2m for plane-wave basis set is 
taken as 500 eV . The convergence criterion of self con- 
sistent calculations is 10~ 5 eV for total energy values. 
By using the conjugate gradient method, all atomic posi- 
tions and unitcell were optimized until the atomic forces 
were less than 0.05 eV/A. Pressures on the lattice unit 
cell are decreased to values less than 0.5 kB. To prevent 
interactions between the adjacent supercells a minimum 
of 10 A vacuum spacing is kept. 

To correct the energy bands and band gap values ob- 
tained by LDA, frequency-dependent GWo calculations 
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TABLE I: Calculated results for Group IV elements, their binary compounds and Group III-V compounds having honeycomb 
structure. Stable structures are identified as PL or LB standing for the planar and low-buckled geometries, respectively. The 
values of angle between neighboring bonds, 6; buckling parameter, A; nearest neighbor distance, d; 2D hexagonal lattice 
constant, |ai| = | ci2 1 = a; cohesive energy, E c ; minimum value of the energy gap, Eq calculated using LDA and corrected 
by GWo with the symmetry points indicating where minimum (maximum) of conduction (valence) band occurs; calculated 
effective charges on the constituent cation/anion, Z*, a ; Poisson's ratio, v\ and in-plane stiffness, C, are given. Some of the 
structural parameters are described in Fig. [T] 
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GROUP 


IV 












Graphene 


PL 


120.0 


- 


1.42 


2.46 


20.08 


semimetal 


0.0/0.0 


0.16 


335 


Si 


LB 


116.4 


0.44 


2.25 


3.83 


10.32 


semimetal 


0.0/0.0 


0.30 


62 


Ge 


LB 


113.0 


0.64 


2.38 


3.97 


8.30 


semimetal 


0.0/0.0 


0.33 


48 


SiC 


PL 


120.0 


- 


1.77 


3.07 


15.25 


2.52/KM-3.89/KM 


1.53/6.47 


0.29 


166 


GeC 


PL 


120.0 


- 


1.86 


3.22 


13.23 


2.09/KK-3.56/KK 


2.82/5.18 


0.33 


142 


SnGe 


LB 


112.3 


0.73 


2.57 


4.27 


8.30 


0.23/KK-0.41/KK 


3.80/4.20 


0.38 


35 


SiGe 


LB 


114.5 


0.55 


2.31 


3.89 


9.62 


0.02/KK-0.00/KK 


3.66/4.34 


0.32 


57 


SnSi 


LB 


113.3 


0.67 


2.52 


4.21 


8.72 


0.23/KK-0.69/KK 


3.89/4.11 


0.37 


40 


SnC 


PL 


120.0 


- 


2.05 


3.55 


11.63 


1.18/rK-6.01/rK 


2.85/5.15 


0.41 


98 










III-V 


GROUP 












BN 


PL 


120.0 




1.45 


2.51 


17.65 


4.61/KK-6.57/rK 


0.85/7.15 


0.21 


267 


A1N 


PL 


120.0 




1.79 


3.09 


14.30 


3.08/rM-5.74/rM 


0.73/7.27 


0.46 


116 


GaN 


PL 


120.0 


_ 


1.85 


3.20 


12.74 


2.27/rK-4.58/rK 


1.70/6.30 


0.48 


110 


InN 


PL 


120.0 




2.06 


3.57 


10.93 


0.62/TK-5.35/rr 


1.80/6.20 


0.59 


67 


InP 


LB 


115.8 


0.51 


2.46 


4.17 


8.37 


1.18/LK-2.68/rK 


2.36/5.64 


0.43 


39 


InAs 


LB 


114.1 


0.62 


2.55 


4.28 


7.85 


o.86/rr-i.92/rr 


2.47/5.53 


0.43 


33 


InSb 


LB 


113.2 


0.73 


2.74 


4.57 


7.11 


o.68/rr-i.7i/rr 


2.70/5.30 


0.43 


27 


GaAs 


LB 


114.7 


0.55 


2.36 


3.97 


8.48 


1.29/rK-2.77/rK 


2.47/5.53 


0.35 


48 


BP 


PL 


120.0 




1.83 


3.18 


13.26 


0.82/KK-1.67/KK 


2.49/5.51 


0.28 


135 


BAs 


PL 


120.0 




1.93 


3.35 


11.02 


0.71/KK-1.14/KK 


2.82/5.18 


0.29 


119 


GaP 


LB 


116.6 


0.40 


2.25 


3.84 


8.49 


1.92/TK-3.50/KM 


2.32/5.68 


0.35 


59 


AlSb 


LB 


114.8 


0.60 


2.57 


4.33 


8.04 


1.49/KM-2.01/KK 


1.58/6.42 


0.37 


35 


BSb 


PL 


120.0 




2.12 


3.68 


10.27 


0.39/KK-0.23/KK 


3.39/4.61 


0.34 


91 



are carried out.— Screened Coulomb potential, W, is kept 
fixed to initial DFT value Wo and Green's function, G, is 
iterated five times. Various tests are performed regarding 
vacuum level, kinetic energy cut-off potential, number of 
bands, k-points and grid points. Final results of GWo 
corrections are obtained by using (12x12x1) k-points in 
BZ, 15 A vacuum spacing, default cut-off potential for 
GWo, 160 bands and 64 grid points. 

Cohesive energies (E c ) per pair of atoms (see Table [TJ) 
are calculated by using the expression 

E c = E T [AB] - E T [A] -E T [B] (1) 

where Ej~ [AB] is the total energy per A-B pair of the 
optimized honeycomb structure; -Et[A] and Et[B] are 
the total energies of free A and B atoms. All of them 



are calculated in the same cell. For graphene, Si and 
Ge, A = B. For the charge transfer analysis, the effec- 
tive charge on atoms are obtained by Bader method^ 
In fact, various methods for charge transfer analysis give 
similar trends for the honeycomb structures studied in 
this paper. For rigorous test of the stability of fully re- 
laxed honeycomb structures under study, we also calcu- 
lated phonon modes by using force constant method^ 
Here the Dynamical Matrix was constructed from forces, 
resulting from displacements of certain atoms in (7x7) 
supercell, calculated by VASP software. 
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FIG. 2: (Color online) Variation of total energy of C, Si, Ge, 
BN, GaAs honeycomb structures with respect to the lattice 
constant a of 2D hexagonal lattice. The stable local minima 
of the Born-Oppenheimer surface for each structure is shown 
with a dashed line separately. 



III. ATOMIC STRUCTURE 

We first present a detailed analysis of two-dimensional 
(2D) hexagonal structure of binary compounds of Group 
IV elements, their binary compounds and Group III-V 
compounds all forming honeycomb structure. In our 
study, we also include the discussion of graphene, Si and 
Ge in honeycomb structure for the sake of comparison. 

Graphene has a 2D hexagonal lattice in which C atoms 



are arranged to form a planar (PL) honeycomb struc- 
ture as shown in Fig[TJ Accordingly, it has a six fold 
rotation axis, Cq at the center of the hexagon, which 
is perpendicular to the atomic plane. Hexagonal lattice 
has a two-atom basis in the primitive unit cell, corre- 
sponding to A- and B-sublattices. That is three alter- 
nating atoms of each hexagon belong to one of the two 
sublattices. In graphene planar geometry is assured by 
the formation of strong 7r-bonding between two nearest 
neighbor p z -orbitals perpendicular to the graphene plane. 
The resulting n- and 7r*-bands determine also relevant 
electronic properties. In addition, there are strong yet 
flexible, covalent cr-bonds derived from the planar hybrid 
sp 2 orbitals between adjacent C atoms. Nearest C atoms 
are separated by 1.42 A and the magnitude of the hexag- 
onal Bravais lattice vector is 2.46 A. Briefly, the planar 
sp 2 hybridization and perpendicular p z orbitals underlie 
planar geometry, unusual mechanical strength and elec- 
tronic structure of graphene. 

In Fig. [2] we show the variation of the total energy 
with respect to the lattice constant a of the 2D hexago- 
nal Bravais lattice. We see that C and BN stayed planar 
and have a single minimum. The situation with Si, Ge 
and GaAs is different, since they have two other min- 
ima corresponding to low buckled (LB) and high buck- 
led (HB) geometries in addition to planar geometry. In 
fact, the total energies corresponding to the minimum of 
the planar geometry are already higher than those of LB 
and HB geometries. In buckled geometries, while atoms 
of A-sublattice are rising, those of B-sublattice are low- 
ered. At the end atoms of A- and B-sublattices lie in 
different planes having a buckling distance, A as shown 
in Fig. [TJ The value of A in HB geometry is high and is 
in the ranges of ~ 2. 5 A, but it is low in LB geometry and 
ranges between 0.4 and 0.7A. We note that two minima 
corresponding to HB and LB geometries in Fig. [2] are 
seperated by a significant energy barrier. As we discuss 
in subsection III.B the minimum of HB is not actually a 
local minimum on Born-Oppenheimer surface. The six- 
fold rotation symmetry of graphene is broken as a result 
of buckling and changes to the three-fold rotation sym- 
metry C3. The similar symmetry breaking takes place 
also in Group III-V compounds having PL honeycomb 
structures. 

In concluding this discussion, we point out that even 
if the calculated total energy has a minimum relative to 
a specific structural parameter, this may not correspond 
to a local minimum. Then it remains to answer which 
of these minima in Fig. [5] corresponds to a local minima 
on the Born-Oppenheimer surface. At this point, reliable 
tests for stability of structure have to be performed. 



A. Phonon modes and stability 

Analysis of phonon modes provides a reliable test for a 
structure optimized conjugate gradient method. If there 
is an instability related with a phonon mode with k in 
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BZ, the square of frequency, fi(k) obtained from the dy- 
namical matrix becomes negative implying an imaginary 
frequency. Phonon calculations are performed by taking 
into account the interactions in (7x7x1) large supercells 
consisting of 98 atoms. For all the infinite 2D honey- 
comb structures there are three acoustical (A) and three 
optical (O) modes. In Fig. [3J we present the calculated 
dispersions of phonon modes of various honeycomb struc- 
tures; namely Group IV elements and their binary com- 
pounds, as well as Group III-V compounds. According 
to these results some materials, for example Sn cannot be 
stable in the honeycomb structure. None of the honey- 
comb structures are stable in the HB geometry. Only 22 
honeycomb structures studied in this paper display £!(k) 
dispersions, which confirm their stability either in PL or 
in LB geometry. 

Well known linear behavior of phonon dispersions of 
LA and TA branches and quadratic behavior of ZA 
branch around T-point also exist for all 2D lattices. 
Among these, the LA and TA phonon branches are heat 
carrying modes. However, it was shown that bend- 
ing branch ZA makes negligible contribution to thermal 
conductivity— Note that only the ZA branch gets imag- 
inary frequencies for k wave vectors in the vicinity of T- 
point in BZ. While this indicates the structural instabil- 
ity for infinitely large periodic structures, it can be taken 
also as evidence that an infinitely large 2D hexagonal lat- 
tice can be stabilized by defects or ripples with large wave 
length. Our analysis for such systems treated in (n x n) 
supercells showed the existence of ripples. Clearly, finite 
patches of such structures can become stable. This argu- 
ment has been confirmed by our calculation of vibration 
modes of finite size patches. 

In our earlier paper— we found that ZA mode of 2D 
periodic Ge honeycomb gets imaginary frequencies near 
T-point of BZ. This situation has been interpreted as 
the instability against long wavelength transversal waves, 
which can be stabilized by rippling or by limiting the size 
of Ge sheets. An extensive analysis of phonon modes in 
the present study revealed that the extent of the region 
of imaginary frequencies around the T-point also depends 
on the mesh size used in the calculations. Decreasing the 
mesh size may lead to the decreasing of their particular 
zone. Therefore, a tedious analysis of the right mesh size 
is required to determine whether or not the imaginary 
frequency zone of ZA mode is an artifact of the mesh 
size. We performed this analysis of mesh size for Ge-LB 
structure and found an optimum fine mesh size where 
imaginary frequencies of ZA mode disappeared. We also 
note that since the interatomic forces related with ZA 
modes decay rapidly, the numerical inaccuracy in calcu- 
lating forces due to the transversal displacement of dis- 
tant atoms may give rise to difficulties in the treatment 
of ZA modes. Briefly, caution has to be exercised in de- 
ciding whether the imaginary frequencies of ZA modes 
is an artifact of numerical calculations. Of course, if the 
presence of imaginary frequencies near the T-point is a 
reality, the instability of infinite and 2D periodic struc- 



tures can be stabilized by finite size or large wavelength 
ripplings. 

Calculated phonon dispersion of graphene is in good 
agreement with previous LDA results and also with re- 
ported experimental dataj 56 i 57 ' 58 Around 1600 cm -1 , LO 
and TO eigenmodes are degenerate at T point. In-plane 
TA and LA eigenmodes have linear dispersion around the 
r point. As it is mentioned in earlier works on 2D struc- 
tures, out-of-plane ZA eigenmode have quadratic phonon 
dispersion in the vicinity of T-point. Here, the calculated 
value for out-of-plane optical eigenmode ZO is around 
900 cm -1 . Existence of strong electron-phonon coupling 
in TO eigenmode at the if-point and E 2g modes at T 
-point is the reason of the Kohn anomaly at these points. 
Therefore, scattering by phonon with the energies that 
corresponds to these modes can cause noticeable decrease 
in transmission spectrum i 58 ' 59 Since force constant de- 
creases with increasing atomic number or row number in 
the Periodic Table, calculated vibration frequencies ex- 
hibit the same trend. 

As a result of symmetry in honeycomb structures of 
Group IV elements (such as graphene, Si and Ge), ZO 
and TO branches cross at K point. We also note that the 
ZO branch of a binary compound comprising at least one 
element from the first row falls in the frequency range of 
acoustical vibration modes. By comparing the phonon 
dispersions of InN, InP, In As and InSb samples, it is 
seen that ZO mode have increasing tendency to move 
apart from the LA and TA modes with increasing nearest 
neighbor distance. However, in all the samples contain- 
ing first row elements, ZO mode is located between the 
LA and TA modes. These characteristic trends of ZO 
mode exists for all the 2D honeycomb structures. 



B. Stability via Puckering 

According to analysis of stability based on the cal- 
culated phonon modes, structures which do not contain 
first row elements occur in LB (puckered) structure corre- 
sponding to a local minimum in Born-Oppenheimer sur- 
face. Through puckering the character of the bonding 
changes. Different hybrid orbitals underlie the different 
allotrophic forms of C atom. While the bonding of dia- 
mond structure is achieved by tetrahedrally coordinated, 
directional sp 3 hybrid orbitals, sp 2 +p z and sp + p x +p y 
hybrid orbitals make the bonding in graphene and cu- 
mulene (monatomic chain of carbon atoms), respectively. 
In forming hybrid orbitals one of two valence s states is 
excited to p state, whereby a promotion energy is imple- 
mented to the system. However, by s and p hybridiza- 
tion the hybrid orbitals yield the maximum overlap be- 
tween adjacent C-C atoms and hence the strongest pos- 
sible bonding. This way, the promotion energy is com- 
pensated and the system attains cohesion. In sp 3 hybrid 
combination one s orbital is combined with p x , p y and 
p z orbitals to form four orbitals directed from the central 
C atom towards its four nearest neighbors in tetrahedral 



6 



Oraphcnc 



Si-LB 



Ge-LB 



BINARY COMPOUNDS a 




FIG. 3: Calculated vibration frequencies of phonon modes SI versus k of Group IV elements, their binary compounds and 
Group III-V compounds having honeycomb structure. Compounds having at least one constituent from first row elements 
have tendency to form planar structure. PL and LB stand for the planar and low-buckled (puckered) honeycomb structures, 
respectively. 
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directions. The angle between these bonds is ~ 109.5°. 
In sp 2 one s is hybridized with p x and p y orbitals to make 
three planar sp 2 which are directed from the central C 
atom at the corners of the hexagons to its three near- 
est neighbors. For the cumulene s orbital is hybridized 
with p z orbital along the chain axis. In this respect, 
the strengths (i.e. self-energy) of these hybrid orbitals 
decreases with increasing number of p-type orbitals in 
the combination; namely sp is strongest whereas sp is 
least strongest. As for the the distance of C-C bonds, 
it is shortest in cumulene (1.29), but longest in diamond 
(1.53 A). In addition to these hybrid orbitals, dangling p 
orbitals make also 7r-bonding between two C atoms. The 
7r-bonding between two adjacent C atoms in graphene 
and cumulene assures the planarity and linearity, respec- 
tively. In Figf3]the charge density of the 7r-bonds between 
neighboring C atoms explains how the stable planar ge- 
ometry is maintained. 

As the bond distance between two nearest neighbor 
atoms increases the overlap of the p z orbitals decreases. 
This, in turn, decreases the strength of the 7r-bond. This 
is the situation in the honeycomb structures of Si, where 
the Si-Si bond distance (2.34 A) increased by 92% rela- 
tive to that of the C-C bond. As a result of weaker it 
bonds in Si can not maintain the stability of the planar 
geometry, and the structure attains the stability through 
puckering, where three alternating atoms of a hexagon 
raises as the remaining three is lowered. At the end 
the structure is buckled by A. Through buckling the 
sp 2 hybrid orbital is dehybridized and s, p x , p y orbitals 
are then combined with p z to form sp 3 hybrid orbitals. 
While three sp 3 hybrid orbitals form covalent bonds with 
three nearest neighbor atoms, one sp orbital directed 
upwards perpendicular to the atomic plane form a weak 
bond with the adjacent sp 3 orbital directed downwards. 
Eventually, a puckered Si honeycomb structure is remi- 
niscent of graphan o 60 i 61 ' 62 with alternating C atoms sat- 
urated with hydrogen atoms from different sites. Briefly, 
puckering occurs as a result of weakening of 7r-bonds, 
whereby the structure regains its stability through tetra- 
hedrally coordinated sp 3 -like bonding. Puckering may 
be explained in terms of Jahn- Teller theorem^ predict- 
ing that an unequal population of degenerate orbitals in 
a molecule leads to a geometric distortion. This way the 
degeneracy is removed and the total energy is lowered. 

As shown in Table []] 11 out of 22 honeycomb struc- 
tures prefer planar geometry, the rest is puckered to re- 
gain stability. There are interesting examples for pla- 
nar and puckered ring structures: Besides BN, planar 
B(CHs)3 molecule and well known B(OH)^, boric acid 
planar crystals are other examples for boron including 
materials. S n rings such as Sr, Sg and S±2 have also 
puckered structures with a crown shaped Another ex- 
ample is planar cyclobutadiene (C4H4); the well known 
molecule was shown that it lowers its energy when it 
has puckered (butterfly) structure for its positive dian- 
ion (CaH^ 2 ). Additionally, it was shown before^ that 
puckered shape of cyclooctatetraene CsHg takes planar 



Planar (PL) Low Buckled (Puckered) 




FIG. 4: In the honeycomb structures of C, Bn and several 
others, the planar geometry is maintained by the strong n- 
bonding through the perpendicular p z orbitals, in addition to 
the a bonding through the sp 2 hybrid orbitals. In the case 
of honeycomb structures formed by the elements beyond the 
first row, the tt-ti bonding is weakened due to the increasing 
bond length. The structure is stabilized by puckering where 
sp 2 hybrid orbitals are slightly dehybridized to form sp 3 -like 
orbitals. This situation is depicted for Si honeycomb struc- 
ture. Increasing charge density is plotted with colors from 
yellow (light) to red (dark). 

shape for its negative dianion CsHg 2 . For the phosp- 
hazenes (NPX2), while N4P4FS is planar N4P4CIS and 
N^P^CHs)^ have buckled shape o 66 ! 67 and thus the rule 
regarding the compounds including first row elements is 
still valid. 

Having discussed the general aspects, we now concen- 
trate on the optimized atomic structure and correspond- 
ing electronic properties of stable honeycomb structures. 
Calculated values of atomic and electronic structural pa- 
rameters are given In Table |U One notes that 11 struc- 
ture having planar geometry has at least one constituent 
from the first row elements of the Periodic Table; namely 
C, B and N. Since the radii of these atoms are relatively 
small, their presence as one of the constituents assures 
that the bond length is small enough to keep strong tt- 
bonding. This explains how the radius of constituent 
atoms enters as a crucial ingredient in the structure. The 
rest of the honeycomb structures in Table |T] including Si 
and Ge are puckered to have LB geometry. 

Finally, we note that calculated results given in Table 
U display interesting trends depending on the radius of 
constituent elements or their row number in the Periodic 
Table. For example, the bond strength or cohesive energy 
E c of a honeycomb structure gets weaker as the atomic 
radii or the row number of the constituent elements in- 
crease. Also band gap Eq and lattice constant a show 
similar trends. 



IV. MECHANICAL PROPERTIES 

Honeycomb structure with sp 2 bonding underlies the 
unusual mechanical properties providing very high in- 
plane strength, but transversal flexibility. We note that 
graphene and its rolled up forms, carbon nanotubes are 
among the strongest and stiffest materials yet discov- 
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ered in terms of tensile strength and elastic modulus. 
We investigated the mechanical properties of 22 stable 
honeycomb structures listed in Table HI We focused on 
the harmonic range of the elastic deformation, where the 
structure responded to strain e linearly. We pulled the 
rectangular unit cell in x- and y-directions in various 
amounts and generated a mesh of data corresponding 
to the strains in x and y directions versus strain energy 
defined as E s — £?r(e) — Et{c — 0); namely, the total 
energy at a given strain e minus the total energy at zero 
strain. The data is fitted to a two-dimensional quadratic 
polynomial expressed by 

Es(e x ,e y ) = aiel + a 2 el + a 3 e x e y (2) 

where e x and e y are the small strains along x- and y- 
directions in the harmonic region. Owing to the isotropy 
of the honeycomb structure a\ = ai. The same equation 
can be obtained from elastic tensor— in terms of elastic 
stiffness constants, namely a\ — ai — (h ■ Aq/2) ■ C\\ ; 
03 = Qi-Aq) ■ C\2- Hence one obtains Poisson's ratio v = 
-Hrans/taxiah which is equal to C 12 /Cn = a z /2a x . Simi- 
larly, the in-plane stiffness, C — h- C\\ ■ (1 — (C11/C12) 2 ) 
= (2ai-(a3) 2 /2ai)/(Ao). Here h and Aq are the effec- 
tive thickness and equilibrium area of the system, respec- 
tively. In Table [H calculated Poisson's ratio and in-plane 
stiffness results are shown. The calculated value of the 
in-plane stiffness of graphene is in agreement with the ex- 
perimental value of 340 (N/m)i§£ Graphene has highest 
in-plane stiffness and lowest Poisson's ration among all 
honeycomb structures of Group IV elements and Group 
III-V compounds. Being a compound of first row ele- 
ments, BN has second highest C and second lowest v. 
The Poisson's ratio v increases with increasing row num- 
ber of elements of elemental and compound honeycomb 
structures. C shows a reverse trend. The order of values 
of C in the last column of Table Q] is similar to that of 
cohesive energies E c in the seventh column. This clearly 
indicates a correlation between E c and C as shown in 
Fig. [5] 

V. ELECTRONIC STRUCTURE 

Our results on the electronic band structure of Group 
IV elements and binary compounds between different 
Group IV elements and Group III-V elements, which are 
stable in either infinite periodic form or in finite size, are 
presented in Fig. [5] In these hexagonal lattice structures 
(PL or LB) relevant electronic energy bands around the 
Fermi level are derived from n and 7r*-bands. In elemen- 
tal honeycomb structures, such as graphene, Si and Ge, 
these bands have linear crossings at two in-equivalent K- 
and i-T'-points of BZ, called Dirac points and hence they 
are semimetallic. Because of their linear dispersion of 
E(\l), the charge carriers near the Dirac points behave 
as massless Dirac Fermions. By fitting the 7r-and tt*- 
bands at k = K + q to the expression, 




10 15 20 

Cohesive Energy (eV) 



FIG. 5: (Color online) A plot showing the correlation be- 
tween the cohesive energy E c and in-plane stiffness C among 
honeycomb structures. 



£(q) ~ v F h|q| + 0(q 2 ) (3) 

and neglecting the second order terms with respect to 
q 2 , one can estimate the Fermi velocity for both Si and 
Ge as vp ~ 10 6 . We note that vf calculated for 2D 
LB honeycomb structures of Si and Ge are rather high 
and close to that calculated for graphene using the tight- 
binding bands. It is also worth noting that because of 
the electron-hole symmetry at K- and i^'-points of BZ, 
2D LB Si and Ge are ambipolar for E(q)= E F ± SE, 8E 
being small. 

In graphene Dirac fermions have a high Fermi ve- 
locity, Vf — c/300. Due to its high carrier mobility, 
graphene based ballistic transistors operating at room 
temperature have already been fabricated^ In addi- 
tion to these unusual electronic properties of graphene, 
the observation of anomalous quantum Hall effect and 
the possibility of Klein paradox are features, which at- 
tract the interest of researchers. Electronic properties 
of graphene and graphene-based structures have recently 
been reviewed] 71 1 72 

In the polar structures, such as BN, GaAs, after charge 
transfer p z orbital electrons are located predominantly on 
one type of atom. Thus the degeneracy of valence and 
conduction bands at K point is removed and gap opening 
occurs* 7 ^ In Table U the minimum width of band gaps cal- 
culated with LDA are given together with the symmetry 
points where the maximum (minimum) of valence (con- 
duction) bands occur. Values of these band gaps after a 
correction by the GWo method are also given. The bands 
of compounds before and after GWo correction are also 
illustrated in Fig. O 
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FIG. 6: (Color online) Energy bands of monolayer honeycomb structures of Group IV elements and binary compounds between 
different Group IV elements and Group III-V elements. All the Group IV elements have semimetallic electronic structure. 
Band structures show that like graphene n and it* bands of Si and Ge have linear band crossing at the Fermi level. Binary 
compounds of Group IV and Group III-V elements are semiconductors. Corrections LDA band gaps using GWo are indicated 
by small circles. Band gaps are shaded. 
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Binary compounds have polar character in addition 
to the covalency of bonds. Effective charge on cation 
and anion Z*^ a , charge transferred from cation to anion, 
5p = Z* — Z v (Z v being the valency of the constituent 
atom) are calculated using Bader analysis. In spite of the 
ambiguities in finding the true effective charge, the cal- 
culated effective charges in Table U give some idea about 
the direction of charge transfer and ionicity of the hon- 
eycomb structure. For some binary compounds like SiC, 
BN, AIN calculated effective charges appear to be right 
in sign but exaggerated in magnitude. We note that as 
the difference in the row numbers of constituent elements 
increases, Z* usually decreases. One can also generalize 
that the charge transfer decreases with increasing row 
number or atomic radii of anion if the cation is fixed. 
This trend is obvious in the structures of InN, InP, In As 
and InSb. 



VI. HETEROSTRUCTURES 

Depending on the constituent elements the band gaps 
of compound honeycomb structures change in a wide en- 
ergy range. In contrast, the lattice constant a of the 
compounds do not show significant variation. The situa- 
tion, where band gaps of two honeycomb structures are 
significantly different, while their lattice constants are 
practically the same, is a convenient condition to make 
semiconductor heterostrucures. As an example, let us 
consider AIN and GaN, which have LDA band gaps of 
3.08 eV and 2.27 eV, respectively. Their lattice con- 
stants are not significantly different and are 3.09 and 3.20 
A, respectively. Moreover, armchair nanoribbons can 
form pseudomorphic heterostructure with perfect junc- 
tion. This is reminiscent of an AIN/GaN commensurate 
heterostructure having 2D interface. Owing to charge 
transfer between constituent nanoribbons at the junc- 
tion, the bands are shifted and eventually aligned. Het- 
erostructures of elemental and compound semiconduc- 
tors generating a 2D electron gas and devices produced 
therefrom have been an active field of study in device 
physics in the past decades. It is expected that the het- 
erostructure of armchair nanoribbons of GaN and AIN 
can constitute a ID analog. When periodically repeated, 
this heterostructure can form superlattices behaving as 
multiple quantum wells or quantum dots. Earlier similar 
effects have been investigated for the heterostructures of 
graphene nanoribbons with different widths^ 

Superlattices of armchair honeycomb nanoribbon 
structures can be constructed according to the width and 
repeat periodicity of the constituent segments. We can 
label GaN/AIN superlattices as GaN/AlN(m,n2;si,S2)- 
Here, Si and S2 specify the length of segments in terms of 
the numbers of the unitcells of constituent nanoribbons. 
Also n\ and n-i specify the width in terms of the num- 
ber of dimer lines in the primitive unit cell of constituent 
nanoribbons. By varying the n and s , we can construct 
variety of superlattice structures. As a proof of concept, 




(1.1) BANDS (C) CHARGES 




FIG. 7: (Color online) (a) A superlattice formed by peri- 
odically repeating heterostructure of armchair nanoribbons 
of GaN-AIN compounds. Red, yellow, blue and small balls 
represent Ga, Al, N and H atoms, respectively, (b) Energy 
band structures of constituent GaN, AIN nanoribbons and 
resulting superlattice in momentum space. The band gaps 
are shaded by yellow, (c) Band decomposed isosurface charge 
densities for lowest two conduction band and highest two va- 
lence bands. Charges of lowest conduction band and highest 
valence band states are confined in GaN side of the junction, 
which has relatively smaller band gap than AIN. 



we consider a superlattice GaNAlN(10,10;4,4) as shown 
in Fig. [7] In the same figure we also presented the 
electronic band structure of constituent GaN and AIN 
nanoribbons. Upon construction, the atomic structure 
is fully optimized. Resulting energy band structure and 
charge density isosurfaces are presented in the same fig- 
ure. 

The highest valence band and the lowest conduction 
band states are flat and they are identified as confined 
states. As a result, one can deduce a type-I (normal) 
band alignment since states are confined to the GaN part 
of heterostructure. One notes that the bandgap of the 
superlattice in momentum space is different from those 



11 



of constituent nanoribbons and can also expect that the 
superlattice band gap in momentum space gets larger as 
the extension of GaN and A1N sides increases. Similar 
confined states can be obtained by constructing A1N core 
and GaN shell structures, where electrons are expected 
to be confined in the core region. 

VII. DISCUSSIONS AND CONCLUSIONS 

In view of the exceptional electronic, magnetic and 
mechanical properties of recently synthesized graphene, 
questions have been raised whether well-known materials 
in micro and optoelectronic industry can attain similar 
honeycomb structures. It is hoped that unusual proper- 
ties can be attained from these structures. The present 
paper examined a large number of materials, Group IV 
elements, binary compounds of these elements, as well 
as a large number of Group III-V compounds to reveal 
whether they may form 2D honeycomb structure. For 
several decades the bulk crystals of these materials have 
dominated micro and optoelectronic industry. Based 
on ab-initio structure optimization and calculations of 
phonon modes we are able to determine 22 honeycomb 
structure, which can be stable in a local minimum on the 
Born-Oppenheimer surface as either 2D infinite periodic 
crystals or finite size flakes (patches). Our calculations 
reveal that Group IV elements, Si and Ge, and binary 
compounds SiC, GeC, SnC, SnSi, SnGe and SiGe have 
stable honeycomb structures. However, while SiC, GeC, 
SnC are planar like graphene and BN, Si, Ge, SnSi, SnGe 
and SiGe are buckled (or puckered) for stabilization. We 
also find that all III-V compounds containing first row el- 
ements B, C or N have planar stable structures. However, 
the binary compounds formed from the combination of 
Al, Ga, In and P, As, Sb are found to be stable in low 
buckled structure. 



For honeycomb structures which were deduced to sta- 
ble, an extensive analysis have been carried out to deter- 
mine their atomic structure, elastic and electronic prop- 
erties. While Si and Ge are semimetallic and linear 
band crossing at the Fermi level like graphene, all the 
binary compounds are found to be semiconductors. In- 
terestingly, these honeycomb materials exhibit interest- 
ing trends regarding cohesive energy, band gap, effective 
charge, in-plane stiffness, Poisson's ratio depending on 
the row numbers of their constituent elements or their 
radius. Interestingly, we deduced a relation between in- 
plane stiffness and cohesive energy among all honeycomb 
structures studied in this work. 

These materials in honeycomb structure have a variety 
of band gaps. Even more remarkable is that the nanorib- 
bon forms of these materials provide diverse properties 
depending on not only their constituents, but also their 
chirality and width. All these properties are expected 
to offer number of applications. Therefore, the studies 
related with their functionalization by vacancy defects 
or adatoms, their mechanical and spintronic properties, 
their heterostructures and core-shell structures will open 
a new field of research. We hope that the findings in this 
work will promote the research aiming at the synthesis 
of these materials. 
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